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Abstract 



The non-linear relation between electric polarization and electric field strength is achieved through intro- 
ducing the retarded electromagnetic interactions between classical charge particles. The result agrees with 
the phenomenological theory in current non-linear optics, means that the non-linear phenomena in optics 
come from the retarded electromagnetic interaction between charged particles. 

The current non-linear optics describes interaction between light and materials by the half-classical 
method, i.e., use the classical theory of electromagnetic field to describe light's motion and quantum theory 
to describe interaction between photons and material particles. The theory can deal with non-linear optics 
well. The Maxwell equation of electromagnetic interaction in medium is 

V*S=f VxH = 7 / + f (!) 

V-D = p f V-S = (2) 

The experiments show that for general non-iron electric medium and non-iron magnetic medium, under the 
condition of weak field, the relations between electric polarization P and electric field strength, as well as 
magnetic polarization M and magnetic field strength are linear, connected by the so-called constructive 
equations 

D = eE = e E + P B = fiH = fio(H + M) (3) 

P = eoXeE M = Xm H (4) 

But the experiments also show that in the strong fields, the relations between P and E, as well as M and 
H are non-linear and vary complex. This kind of relations can't be deduced from theory at present. We 
now get them by phenomenological hypothesis. For example, in the non-linear optics, the relation between 
electric polarization P and electric field strength is supposed to be W 

P = e (xlE + x 2 e -EE + £:EEE--) (5) 

Here x\ ' IS polarization tensors. By means of the formula and the Maxwell equation, so much non-linear 
phenomena in optics can be described well. It is proved blow that after the retarded electromagnetic 
interaction between charged particles are considered, we can obtain Eq.(5) easy and shows that the non- 
linear phenomena in optics come from the retarded electromagnetic interaction between charged particles. 

Let t'jfjV' and a' represent retarded time, coordinate, velocity and acceleration, t,r,v and a represent 
non-retarded time, coordinate, velocity and acceleration. A particle with charge % velocity Vj and acceler- 
ation a'j at space point Tj(t') and time f would cause retarded potentials as follows at space point r.i(t) and 
time t 

v > = %^ — = — — (6) 

4**0(1 4**0(1 
In the formula fy(t,f ) — rj(t) — rj(t'), =| |, = f«/rL. When the particle's speed v « c, the 
retarded distance r^-(t') at retarded time t' can be replaced approximately by non-retarded distance r.y(i), 
i.e., we can let 

t' =t-r' lJ (t')/c^t-r l0 {t)/c (7) 
rij(t, =| n{t) - i$(f) 1 = 1 n(t) - - 4-(f )/c] H| n{t) - r*-[t - rij (t)/c] |= r'^t) (8) 



It is noted that r£,-(t) i= fij(t), for r^-(i) is the approximate retarded distance, but rij(t) is not the retarded 
distance. In this case, we can develop retarded quantities into series in light of small quantity rij/c. By 
relation vj = dfj/dt, we get ^ 

r--(M) - n(t) r-(t - nj /c) = n(t) fj(t) + ^r«(t) - ^>r&(t) + ^^(t) + ■■■ 



= -. (/) + M rtf(t) _M r&(t) + M r 3 ((t) + ... 



er^f ) * ^(t) - (t) + + • • • aj(f ) ~ a^t) - ^r«(t) + • • • (10) 

Mediums are composed of atoms and molecules, and atoms and molecules can be regarded as electrical 
dipolcs composed of two particles with charges ±g. Suppose the distance between two charges is (t) at time 
t, the electrical dipole moment is Pj. The direction of Pj = Qjfij is from the point the j -particle located 
pointing to the mass center point fj of electric dipole. By considering the retarded effect of interaction, 
according to Eq.(9), the retarded electrical dipole moment becomes approximately 

P'j = Qi^a = Qjinj + + ^Lr%) (11) 

Vj(t)3j(t) and Sj(t) are the -particle's velocity, acceleration and acceleration of acceleration individually. If 
there is only one electrical dipole in space, the directions of fy, Vj, ctj and Sj are on the same straight line. 
In this case, the subscript can be removed and Eq.(ll) can be written as 

i y -(Kr + 2r-^r» + A r 3 + ... ) (12) 

When there exists the action of external electric field E = E sm(uit + 0' ), the motion equation of electrical 
dipole is 

The second item on the right side is damping force (radiation damping force and damping force caused by 
collisions between particles). When w 2 > 0, the solution of Eq.(13) is 



r = Asm{ut - ) + Be~ 0t sm(^J u 2 - pH - 5) (14) 
A= - qE ° e =e' +tg-\^-^) (15) 

When t is big enough, we have e^* 3 '/ 2 — > and the second item on the right side of Eq.(14) can be omitted. 
In this case, we have 

v = Alo cos(wt — O ) a — —Ato 2 sin(u>i — 60) a = —Aoj 3 cos(wi — 9o) (16) 

Put them into Eq.(12), we get 

2 

P' = qA{sm(cot - 9 Q ) + — sin(wt - 9 ) cos(u;t - 9 ) + ^sin 3 (wt - 9q) 
— — -^sin (ut — 80) cos(cjt — 80) + ■ ■ •} 

3uJ 2 U) U) 2 UJ 2 

= qA{(l + — ) sm(cot -8 ) + —{l- —) sin(2^ - 29 ) - ^ sm(3ujt - 39 ) 

+ ^M^t-46 ) + ---} (17) 

It is obvious that after the retarded interaction is considered, the waves of 2, 3, 4- •• multiple frequencies would 
appear in the vibrations of electrical dipoles. If two electric fields with different frequencies are imported 
with 

E = E 1 sin(wii - 6[) + E 2 sin(w 2 t - 9' 2 ) (18) 



similar results can be obtained 

r = A x sinf^xi - 6{) + A 2 sin(w 2 t - 9 2 ) (19) 

v = Ailui cos(wit — 61) + A2LO2 cos(w2i — 62) (20) 

a = —A\uj\ sin(wii — 6\) — A2UJ2 sin(w2i — #2) (21) 

a = — Aiujf cos(wii — 6*i) — A 2 u>2 cos(w2^ — #2) (22) 



In the formula, A { = gB /m\/( w o - ^fY + ( 2 A^) 2 , ^ = ^ + ^ _1 [ 2 /^/(wo ~ uj2 )}- Put them into E q-(12), 
we get 



P> = qiiA, + + A^K^ + u 2 )] sin( Wl t - 00 + [A 2 + ^ + .4 2 A 2 (c 2 + ^)] sin(c 2 t - 

+ ^2) sin((wi + w 2 )* - (#i + 2 )) + (0J1 - u> 2 ) sin((o) 2 - wi)t - (8 2 - 61))] 

2c 



1) 



H — i-^-[sin(3wit - 36>i) - 3sin(wii - 6>i)] + — ^[sin^i - 36> 2 ) - 3sin(w 2 t - #2)] 



4M 

4c L ~ ~~ K " ±/J ' 4c 

w 2 

-A x Al{-^ + w^)[sin((wi + 2uj 2 )t - (0i + 26> 2 )) + sin((w 1 - 2uj 2 )t - (6>i - 26> 2 ))] 
uj 2 

-A 2 A 2 (uj 2 + -^)[sin((w2 + 2wi)t - (0 2 + 20i)) + sin((cj 2 - 2wi)i - (6 2 - 20i))] + • • •} (23) 

Besides multiple frequencies, sum frequencies and different frequencies would appear. 

Because v = Alo cos(uit — 9 ) = Auj sm(uit — 9 — 7r/2), there is a difference of phase tt/2 between v/u and 
r. There is also a difference of phase tt/2 between a/u 3 and r. If these differences are omitted, we can let 
v ~ ur, a ~ oj 3 r. Considering a = —oj 2 r, put the relations into Eq.(12), we get 

P' ~q(r+-r 2 + ^r 3 + ^r 4 + ---) (24) 
c 2c z 6c d 

According to Eq.(14), when t is large enough and the diffcrcned of phase is omitted, the displacement 
of dipole directs radio to electric field's strength, i.e., r = aE, a = q/m^J (cjq — lo 2 ) 2 + (2(3iv) 2 . For the 
vibrations of atoms in crystal lattice, we have ui » ujo and u >> (3, so a ~ q/mui 2 . Put the relation into 
the formula above, we get 

P^ qa( E + ^E* + ^E> + ^LE> + ... ) ,25) 

This is just the non-linear relation between electrical polarization and electric field strength. Suppose the 
wavelength of incident wave is A = 4 x 10~ 7 M, the dipole vibrating in lattice is hydrogen with mass m, we 
get 

— = — ~ 3.14 x 10- W M/V (26) 
c emu; 

The radio of parameters between the first item and second item is about 10~ 10 . This value coincides with 
experiments ( 4 ). Eq.(27) can be written as 

„/ / 1 „ 9 wet n ut 2 a 2 n u> 3 a 3 / , 
P' c e a (xlE + xl—E 2 + -^E 3 + -^E 4 + ■ ■ ■) (27) 

1 = W 2 = g^« 2 3 = ^C? 4 = g^ 3 ^ 4 (9R] 

eo Soc 2eoc z be^c 6 

But it should be noted that the results exists phase differences between on Eq.(29) and (12) and Eq.(12) is 
fundamental. 

In the practical systems, there exist many dipoles and interactions between dipoles. Besides vibrations, 
dipoles may rotate and shift, so the directions of f{j, Vj, 3j and fij are not on the same straight lines in 



general. By taking sum and let qjdj = Soxtjcn Qj 0J i a 'j/ c = £ oX% a ' "> the total electrical dipole moment of 
medium can be written as the form of partial quantities 

N N 

P-=J2 P ^= £ E^> E + x 2 ej<T E 2 + X % a E 3 + ■■■) (29) 
i=i i=i 

Because the model of dipole is only a kind of approximation, by considering interactions between particles as 
well as by considering other complex factors, in general situation, the non-linear relation between electrical 
polarization and electric field strength should be written as the form of product form of tensors 

P = e a (xlE + x1-EE + £:EEE + ---) (30) 

This is just the foundational relation used in non-linear optics. 

In order to explain the origin of non- linear phenomena in optics or Eq.(32), the non-linear oscillator 
model is discussed at present. The motion equation of oscillator is written as W 

g + 2/4+^r + i> 2 = - (31) 
dt z at m 

In the formula, mDr 2 is the non-linear force to cause non-linear effects. However, it should be seen that the 
origin of this item is still unclear. Where it comes from? Why is takes such form, not mDr 1 ! 2 or mDr 3 
and so on. The current theory can't provide clear explanation. Meanwhile, the value of parameter D can't 
be decided by theory. Quantum mechanics provides a method of perturbation approximation to calculate 
polarizability, but does not yet provide explanation for the origin of non-linear phenomena. So it is proper 
to consider that the origin of non-linear in optics comes from the retarded electromagnetic interaction. 

Because v — > — v, a — > a, a — > a under time reversal, so r 1 shown in Eq.(ll) can't keep unchanged under 
time reversal. In this way, it can be said that most of non-linear phenomenon in optics violate symmetry of 
time reversal. 
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